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Abstract: Modified gravity models such as Hofava-Lifshitz gravity or Einstein-aether theory 
violate local Lorentz invariance and therefore destroy the notion of a universal light cone. 
Despite this, in the infrared limit both models above possess static, spherically symmetric 
solutions with "universal horizons" - hypersurfaces that are causal boundaries between an 
interior region and asymptotic spatial infinity. In other words, there still exist black hole 
solutions. We construct a Smarr formula (the relationship between the total energy of the 
spacetime and the area of the horizon) for such a horizon in Einstein-aether theory. We further 
show that a slightly modified first law of black hole mechanics still holds with the relevant area 
now a cross-section of the universal horizon. We construct new analytic solutions for certain 
Einstein-aether Lagrangians and illustrate how our results work in these exact cases. Our 
results suggest that holography may be extended to these theories despite the very different 
causal structure as long as the universal horizon remains the unique causal boundary when 
matter fields are added. 
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1 Introduction 

In the four decades since the seminal work of Bardeen, Carter, and Hawking [1] on the laws 
of black hole mechanics, a tremendous amount of effort has gone in to understanding hori- 
zon behavior. From the discovery of Hawking radiation [2], and the recognition that the four 
laws have a thermodynamic interpretation [3], to holography [4, 5] and its concrete realization 
through the AdS / CFT correspondence [6-8] , the physics of horizons has provided useful infor- 
mation about quantum gravity. Using horizon thermodynamics and some mild assumptions 
about the behavior of matter, one can reverse the logic and reconstruct general relativity as 
the thermodynamic limit of a more fundamental theory of gravity [9—11]. Integral to horizon 
dynamics is the first law of black hole mechanics, which for the simplest Schwarzschild case, 
and the most similar to what we are interested in, is just 

5M ABM = —— — . (1.1) 

07T(_x N 

Here M ADM is the ADM mass of the spacetime and n Kli and A Kli are the surface gravity and 
cross-sectional area evaluated on the Killing horizon, respectively. Identifying (87tG n ) _1 k kh 
as the temperature of the horizon and the area with the entropy allows one to make the 
analogy with the first law of thermodynamics, 5E = T5S. 
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In order to have an appropriate formulation of thermodynamics for the black hole itself, 
as well as the combined system of the black hole and the exterior environment, it is necessary 
for the horizon to have an inherent entropy. If there was no entropy associated with the 
horizon, one could simply toss objects into the black hole and reduce the total entropy of 
the black hole and exterior system, thereby violating the second law of thermodynamics. If, 
however, the horizon has an associated entropy, then one can reformulate the usual second 
law into the generalized second law 

<5 ("^outside + S'horizon) = 0. (1-2) 

The generalized second law and the thought experiments behind it imply that any causal 
boundary in a gravitational theory should have an entropy associated with it. In general 
relativity, the entropy can be shown to be proportional to the area of a slice of the Killing 
horizon. However, if one includes higher curvature terms the entropy is still a function of 
the metric and matter fields evaluated on a slice of the Killing horizon, though no longer 
proportional to the surface gravity alone [12]. 

A much stronger departure from general relativity comes when one considers models 
that allow vacuum solutions with non-zero tensor fields besides the metric. In these models 
Lorentz symmetry and the equivalence principle are in general broken. A simple example of 
such a model is that of Einstein-aether theory [13], which introduces an aether vector field 
u a and a dynamical constraint which forces u a to be a timelike unit vector everywhere. The 
introduction of the aether vector preserves general covariance, but allows for novel effects 
such as matter fields travelling faster than the speed of light [14] and new gravitational wave 
polarizations that travel at different speeds [15]. Given certain choices of the action for u a , 
the theory can be made phenomenologically viable [16, 17], have positive energy [18], and 
be ghost free [19]. In addition, the aether vector establishes a preferred frame and causality 
can be imposed in that frame [20] by requiring that all matter excitations propagate towards 
the future, even if the momentum vector of an excitation is spacelike. Since there is a 
preferred frame, Lorentz invariance does not hold, nor do the usual Lorentz invariance based 
arguments that a spacelike momentum vector in one frame immediately imply the existence 
of past directed momentum vectors. Thus, the propagation faster than the speed of light does 
not violate causality. 

Even though there is a notion of causality, it seems at first glance as if there would be no 
causal boundaries equivalent to an event horizon in the Einstein-aether theory - by coupling 
the aether vector u a to matter kinetic terms, the matter Lagrangian can be chosen to make 
matter perturbations about flat space propagate arbitrarily fast. This is incorrect - causally 
separated regions of spacetime can exist even in this case. In the static, spherically symmetric 
solutions of Einstein-aether theory found by Eling and Jacobson [21] there exists such a region, 
and the boundary of this region has been dubbed the "universal horizon" 1 . Since the notion 
of a causal boundary and infinite speed modes is counterintuitive, we give a brief explanation 

1 This point was first noted by Sergey Sibiryakov at Peyresq, 2010. 
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of why they can occur here, postponing an in-depth discussion of universal horizons until 
later. 

Consider a static, spherically symmetric spacetime, and cover it with Eddington-Finkelstein 
type coordinates such that the metric takes the form 



and the time translation Killing vector is x a = {1)0,0,0}. Now let denote a surface 
orthogonal to the aether vector u a , so that U is the "aether time" generated by u a that 
specifies each hypersurface in a foliation. At asymptotic spatial infinity \ a an d u a coincide, 
but as one moves in towards r = each Y<u hypersurface bends down to the infinite past in v, 
eventually asymptoting to a 3 dimensional spacelike hypersurface on which (u ■ x) = 0, which 
implies that the Killing vector x a becomes tangent to Y>jj- This hypersurface is the universal 
horizon. It is a causal boundary, as any signal must propagate to the future in U, which is 
necessarily towards decreasing r at the universal horizon. The surface is regular, and in fact 
Barausse, Jacobson and Sotiriou [22] have numerically continued the solution for metric and 
aether fields beyond the universal horizon. 

Since such a causal boundary exists, it is natural to speculate that there must be an 
entropy associated with the universal horizon as well. In spherical symmetry, one does not 
need to worry about the zeroth law of black hole mechanics, as the symmetry enforces that 
all geometric quantities are constant over the universal horizon automatically. Hence one 
can immediately proceed to derive a Smarr formula and corresponding first law. There 
are subtleties, however, as the boundary data of the theory at infinity naively contains two 
parameters. However, as noted by others [21, 22], the boundary data can be reduced to one 
parameter, which is the ADM mass, by requiring that the solution is regular outside the 
universal horizon (see section 4 for details). If one considers only regular solutions, we show 
that there is a Smarr relation between the ADM mass and geometric quantities evaluated 
on the universal horizon. In particular, the resulting Smarr relation contains a contribution 
from the extrinsic curvature of the T<u hypersurface in addition to the standard surface gravity 
term. We further show using a scaling argument that if one considers a transition between two 
regular solutions then there is also a first law that may admit a thermodynamic interpretation. 
With certain choices of the Lagrangian for Einstein-aether theory, we construct new, exact 
solutions and use those as examples to gain insight into the thermodynamics of the first law. 

It is important to note that in previous works, black hole thermodynamics has failed when 
the limiting speed of matter fields is finite, but not necessarily the speed of light. In these cases 
one can construct perpetuum mobiles that violate the second law [23-25]. We will not consider 
any specific matter action in this paper, as our purpose is simply to determine whether a first 
law of mechanics holds for the universal horizon so that a thermodynamic interpretation 
might exist. However it is certainly possible that a true thermodynamic interpretation could 
only hold if the universal horizon is the only causal boundary for all fields. This could, for 
instance, be done if the Lagrangian for any matter fields contained higher derivatives that 
made the local speed of excitations infinite as the energy increased. From the point of view of 




(1.3) 
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effective field theory this is natural as one expects all terms consistent with the symmetries 
of the problem to appear in any operator expansion. 

The paper is organized as follows. We first provide the background for Einstein-asther 
theory in Sec. 2. In Sec. 3 we construct a Smarr formula and first law for spherically sym- 
metric solutions. We then use two new exact, analytic solutions as examples in Sec. 4 to 
verify the Smarr formula and the first law, as well as discuss the regularity of these solu- 
tions. We conclude with some more speculative comments on how to proceed to establish the 
thermodynamic connection and the obstacles that still remain. 

2 The Einstein-aether theory 

Einstein-asther theory was originally constructed [13] as a mechanism for breaking local 
Lorentz symmetry yet retaining as many of the other positive characteristics of general rela- 
tivity as possible. In particular it is the most general action involving the metric and a unit 
timelike vector u a that is two-derivative in fields and generally covariant. General covariance 
is maintained by enforcing the unit constraint on u a via a Lagrange multiplier. Following the 
presentation in [18] the action of Einstein-asther theory is a sum of the usual Einstein-Hilbert 
action 5 EH and the aether action [13] 

5 = 5 EH +5 aE = — -V- f^x^g (2.1) 



WttG 



a- 



In terms of the tensor Z ab cd defined as 2 

Z ab cd = c lg ab g cd + c 2 5 a c 5 b d + c 3 5 a d 5 b c - Ci u a u b g cd , (2.2) 

where Cj, i = 1,...,4 are "coupling constants" (or couplings, for short) of the theory, the 
asther Lagrangian Jzfae is given by 

- J?* = Z ab cd (V a u c )(V b u d ) - X(u 2 + 1). (2.3) 

The aether Lagrangian is therefore the sum of all possible terms for the aether field u a up to 
mass dimension 2 and a Lagrange multiplier A implementing the normalization condition 3 

u 2 = -1. (2.4) 

An additional term, Ht a f,u a u b is a combination of the above terms when integrated by parts, 
and is not included here. 

There exists a number of theoretical as well as observational bounds on the couplings 
Cj, % = 1, . . . , 4 — see e.g. [16, 17] for a comprehensive review. In this work, we assume the 
following constraints to hold on these couplings 

0^ci 3 <l, 0^ci4<2, 2 + ci 3 + 3c 2 >0, (2.5) 
2 Note the indicial symmetry Z ba dc = Z ab cd . 

3 We use a convention where the metric has signature (— , +, +, +)• 
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where we have defined C13 = (ci + C3) and C14 = (ci + C4). As we will see, these combinations 
of couplings, as well as C123 = (ci + ci + C3), play a more direct role in our analysis than 
the individual couplings q. Finally, the constant in the action (2.1) is related to G N , 
Newton's gravitational constant, via 

G m = (l - ^) G N . (2.6) 



2 

as can be established using the weak field/slow- motion limit of the Einstein-asther theory [26]. 
The corresponding equations of motion, obtained by varying the metric, u a , and A, are 

g a b = T l T b , M a = 0, u 2 = -l, (2.7) 

respectively, where the asther stress tensor CT? is given by 

T a f = Xu a u b + c 4 a a a 6 - -gabY c d V c u d + ci [(V a u c )(V b u c ) - (V c u a )(V c u 6 )] + V c X c ab , (2.8) 

and to make the notation more compact we have introduced the following 

M a =V b Y b a +Xu a + c i (V a u b )a b , Y% =Z ac bd V c u d , X\ h =Y\ a u b) - u (a Y b) c + u c Y (ab) . (2.9) 

The acceleration vector a a appearing in the expression for the asther stress tensor is defined 
as the parallel transport of the asther field along the asther field 4 

a a = V u u a . (2.10) 

2.1 Static, spherically symmetric expansions 

We now turn to the static, spherically symmetric case and provide some useful expansions 
which we will later use to analyze the equations of motion. Eling and Jacobson observed that 
for any spherically symmetric solution of Einstein-asther theory, u a is hypersurface orthogo- 
nal 5 [21], which in turn implies that the twist of u a vanishes, i.e. 

u [a V b u c] = 0. (2.11) 

An immmediate consequence of the hypersurface orthogonality condition is that there exists 
a one-parameter redundancy among the couplings Cj. Using the unit norm constraint (2.4), 
the squared twist can be expressed as 

u 2 = (V a u b )(V a u b ) - (V a u b )(V b u a ) + a 2 , (2.12) 

which also vanishes for hypersurface orthogonal solutions. We can, therefore, add any multiple 
of uj 1 to the action without affecting the solutions. In particular, by adding 

A5 = ^ / d 4 x y/^g" w 2 (2.13) 



16vrG a 



4 We use the conventional notation Vx for the directional derivative (X a \7 a ) along any vector field X a . 
Once the normalization condition (2.4) is imposed, the acceleration is always orthogonal to the aether field (i.e. 
u ■ a = 0), and therefore, is always spacelike. 

5 In fact they are also solutions of Hofava gravity [27, 28]. 
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to the action (2.1), where cq is an arbitrary real constant, the sole effect would be to obtain 
a new asther Lagrangian Jzf^,, otherwise identical to (2.3), except with a new set of coupling 
constants i = 1, ...,4, related to the un-primed Cj through 

c[ = ci + co, c 2 = c 2 , c' 3 = c 3 - c , c 4 = c 4 -c . (2-14) 

Thus, by appropriately choosing Co, one can set any one of the couplings ci, C3 and c 4 , or 
any appropriate combinations of them, to any preassigned value. On the other hand, the 
coupling C2, as well as combinations like C13, ct 4 and C123 (2.5) stay invariant under the above 
redefinition of the couplings (2.14). 

Our analysis will be further facilitated by defining a set of basis vectors at every point in 
spacetime so that we can project out various components of the equations of motion. Staticity 
and spherical symmetry implies the existence of a time translation Killing vector x a as we U as 
three rotational Killing vectors i = 1, 2, 3 (only two of the three are linearly independent). 
It is often convenient to choose \ a as one °f the basis vectors, but in this case it is actually 
more helpful to use a different basis. We first take u a to be the (timelike) basis vector. We 
next pick any two spacelike unit vectors, call them m a and n a , both of which are normalized 
to unity, are mutually orthogonal and lie on the tangent plane of the two-spheres 2? that 
foliate the hypersurface Y>jj- Finally, we use s a , the spacelike unit vector orthogonal to u a , 
m a and n a that points "outwards" along a £[/ hypersurface. Note that the acceleration a a 
only has a component along s a by spherical symmetry, i.e. 

a a = (a-s)s a . (2.15) 

Thus, our tetrad consists of {u a , s a , m a , n a }. By spherical symmetry, any physical vector may 
have components along u a and s a , while any rank- two tensor may have components along 
the bi- vectors u a Ub, U( a s b)> u [a s b]> s aSb and g a b, where g a b is the projection tensor onto the 
two-sphere S, bounding a section of a Y>jj hypersurface. For example, the basis-expansion of 
the extrinsic curvature of a T,jj hypersurface is 

K 

K ab = K s a s b + —g ab , (2.16) 

where Kq and K G are scalar parameters related to each other through 

K = K + K, (2.17) 

with K the trace of the extrinsic curvature of the E[/ hypersurface. We ask the reader to 
refer to appendix A for further details on these points. 

6 As a consequence of spherical symmetry and u a being orthogonal on the hypersurface Ey, we have s A ds = 
- compare with (2.11). Thus the vectors s a are orthogonal to the hypersurfaces {E s }, foliating the spacetime. 
It can then be shown that K is the trace of the extrinsic curvature of the two-spheres 23 due to their embedding 
in E s . See appendix A for further details. 
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Finally, we work out some geometric quantities related to the Killing vector x a - We first 
use the spherical symmetry to write the timelike Killing vector \ a as 

X° = -(u-xK + (s-xK . (2.18) 

The basis-expansion of 'VaXb takes the following form 

VaXb = -K(u a s b - s a u b ) . (2.19) 

Here k is defined as surface gravity on any two-sphere 23 since, as we will now show, (at least 
outside the Killing horizon) k is the acceleration of a static observer on "B as measured by an 
observer at asymptotic infinity. First consider the region outside the Killing horizon where 
there exists an outward pointing spacelike unit vector r a , the radial unit vector, orthogonal 
to x a - The tangent vector \ a to the world line of a static observer on any two-sphere 2? is 
simply the unit timelike vector along x a - I n terms of the redshift factor p, we can then write 
X a = PX a (i- e - P 2 = ~X " x)- Using (2.19), the directional derivative of x a along itself is given 
by V^x a = a x r<1 ' w here a x = (n/p) is the local acceleration of the static observer. In other 
words, k = pa x is the redshifted acceleration with respect to an observer at infinity, which 
prompts us to call n the surface gravity on the two-sphere "B. The mathematical analysis 
also follows through inside the Killing horizon 7 and the local acceleration is still given by 
(k/p), but the interpretation of p as the redshift factor no longer holds. Nevertheless, we will 
continue to call k the surface gravity. Using the results from appendix A, the surface gravity 
is explicitly given by 8 



« = ]j-^aXb)(V a X b ) = "(« • s)(u ■ X ) + K (s ■ x) • (2.20) 

As mentioned in the introduction, the universal horizon occurs when the Killing vector 
becomes tangent to a £[/ hypersurface. Therefore as one travels inwards from spatial infinity 
along a T,jj hypersurface, the universal horizon is actually reached only as a limit. Hence 
our quantities defined as "on the universal horizon" refer to this limit, rather than some 
actual intersection of T,jj and the universal horizon which is a hypersurface of constant r 
(again in Eddington-Finkelstein coordinates). On the universal horizon, (u ■ x)uh = an d 
Ox) uh — IxIuh where |x|uh is the magnitude of the Killing vector \ a on the universal horizon. 
Therefore the surface gravity on the universal horizon is 

k uh = K ,uh|x|uh • (2.21) 

7 Note however that inside the Killing horizon \ a is spacelike and r a is timelike. Associating a spacelike 
unit vector with an observer is allowed here since there is no local limiting speed. 

8 We thank Ted Jacobson for pointing out reference [29] in this context, where the surface gravity at the 
Killing horizon is generally shown to be proportional to the expansion of a congruence of timelike geodesies. 
However, we emphasize that in the present paper, equation (2.20) holds everywhere in spacetime rather than 
just at the Killing horizon. We also note that using (u ■ x)kh — —{s ■ x)kh at the Killing horizon, we have 
Kkh = — («■ x)kh{(<i - s) + Ko}kh from (2.20). However, this relation and the the central result of [29], although 
similar in appearance, actually differ since the asther does not define a geodesic flow. 
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By spherical symmetry, the surface gravity is constant over any two-sphere, and thus on the 
universal horizon as well. 

2.2 Equations of motion for the static, spherically symmetric case 

We next study Einstein's equations and the aether equations of motion (2.7) by explicitly using 
the time translational and spherical symmetries of the problem, in addition to hypersurface 
orthogonality. To set up the Einstein's equations, we need to know the basis-expansion of 
the aether stress tensor and the Ricci tensor. By spherical symmetry, they take the following 
form 

lab = ?uuUaU b - 2T*u {a s b) + r^SaSb + ^-g afe , (2.22) 

and 

flab = fiuuUaUb - 2$i us U( a S b ) +2LssS a Sb + ygafe , (2.23) 

respectively. The coefficients of in (2.22) are computed from the general expression (2.8) 
for the stress tensor, using the results in appendix A. The corresponding coefficients for 
%,abi on the other hand, are computed from the defining equation [V a ,V&]X c = — ^. a ^ c ^ d 
by choosing X a = u a or s a , and then contracting the expression again with u a and/or s a 
appropriately 9 . For our present purpose, it is sufficient to show the components T^. and ^ us , 
which are given as follows 

?Z = c u (K(a ■ s) + V u (a • s)), H us = (K - K/2)k - V S K , (2.24) 

where k is the extrinsic curvature of the two-spheres 23 due to their embedding in T,jj (see 
appendix A) 10 . Comparing (2.22) and (2.23), we see that there are altogether four non- 
trivial components of the Einstein's equations. 

The aether's equations of motion (2.7), on the other hand, reduce to the scalar equation 12 

(s-M) = c 13 V s K + c 13 (K Q - K/2)k + c 2 V s K - <T* = , (2.25) 

as a consequence of hypersurface orthogonality and spherical symmetry. Quite naturally, the 
coupling constants that appear in (2.25) above are precisely those which are invariant under 
(2.14). 

A well-known fact in general relativity is that the Bianchi identities (a consequence of 
general covariance) can be used to show that a subset of the Einstein equations are actually 
constraint equations. In Einstein- asther theory, which is also generally covariant, there are 

9 We note that the coefficient cannot be constructed in this method. This is not an obstruction for 
most part since we do not need the explicit expression for In section 4 we use the explicit coordinates to 
construct %^ . 

10 Intimately related to this is the fact that fc/2 is the coefficient of g a i, in the basis-expansion of V a Sf,. 
11 We have one extra Einstein's equation, because we do not assume staticity here. Thus these equations 
can also be used to study time dependent but spherically symmetric perturbations around static solutions. 
12 After solving for the Lagrange multiplier there is no non-trivial projection of M a (2.9) along u a . 
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generalized Bianchi identities, and projections of these give rise to constraint equations as 
well. As explained in [30] (see also [22]), the generalized Bianchi identities for Einstein-asther 
theory are 

V a [g ab - + u a M b ] + M a V b u a = , (2.26) 
and the corresponding constraint equations for the £(/ hypersurfaces are [22, 30] 

($ ab -T < %)u a -M b = 0. (2.27) 

Projecting (2.27) along s b and using (2.24) and (2.25) the explicit form of the constraint 
equation, adapted to the foliation Hy, is 

ci 23 V s if - (1 - c 13 )(K - K/2)k + (1 + c 2 )V s K = . (2.28) 

On the other hand, subtracting the equation Ci3(^. ns — T^j.) = 13 from the aether's 
equation of motion (2.25), we get 

ci 23 V s K = (1 - c n )TZ . (2.29) 

Equations (2.28) and (2.29) are two independent linear combinations of the asther's equation 
of motion (2.25) and the Einstein equation 9^ us = 1™ s . Therefore the two sets of equations 
are equivalent. In section 4, we study these equations along with the uu, ss and the spherical 
components of the Einstein's equations. 

Finally, considering the projection of the aether equations of motion along the Killing 
vector x a > we arrive at the following equation of central importance in this paper (as it will 
be used to derive a Smarr formula) 

V b T ab = 0, T ab = q{u a s h - s a u b ) , (2.30) 

where the quantity q is given by 

q = ~ (l - "y) (a ■ s)(u ■ X ) + (1 - c 13 )K (s ■ x ) + ^K(s ■ X ) . (2.31) 

The derivation of (2.30) closely follows the manipulations leading to the Smarr formula in 
[1]. We also found the algebraic relations of appendix A, especially those discussed in the last 
part of the appendix, useful in arriving at (2.30). 

The similarity between (2.30) and the source free Maxwell's equations allows us to solve 
(2.30) exactly once we adopt a particular coordinate system. A useful choice is Eddington- 
Finklestein like coordinates (see equation (4.1) in section 4), as these coordinates are good 
everywhere in the spacetime. Because of staticity and spherical symmetry, in this particular 
coordinate system T ab has only one non-trivial component, namely T vr . Therefore, solving 
(2.30) amounts to solving for the electrostatic field of a point charge at r = in this particular 

13 This equation, obviously, follows from the Einstein equation {K_ U s — Tus) = 0- 



- 9 - 



geometry. By Gauss's law, we conclude q ~ r 2 . Using the asymptotic series solutions (4.6) 
we then fix the constant of proportionality and obtain 

« = 0-t)&- (2 - 32 > 

where ro is the single free parameter which defines the regular Einstein-asther black hole 
solutions. We ask the reader to refer to section 4 for a more detailed discussion of how the 
parameter counting works and related issues. The parameter ro, as we show in the following 
section, defines the mass of the Einstein-aether black holes. With the aid of equations (2.30), 
(2.31) and (2.32) and the fact that the Einstein- aether black holes constitute a one-parameter 
family of solutions, one can provide very simple derivations of the Smarr relation and the first 
law for Einstein- aether black hole mechanics, which we now turn to. 

3 The Smarr formula and first law 

In this section we give very simple derivations of the Smarr formula and the first law of 
(universal) horizon mechanics for a general static and spherically symmetric Einstein-aether 
black hole. 

To begin with, we need a suitable definition of the mass of the Einstein-asther black holes. 
In this case, the ADM mass of a solution is identical to its Komar mass. From the general 
definition of the Komar mass of a stationary solution 

M ADM = J dE a6 V a X 6 , (3.1) 

23 oo 

where d£ a b is the integration measure on any two-sphere 23, explicitly given by 

dT.ab = -U[ a s b ]dA , 

with dA the differential area element on the two-sphere 23, and 23 oo is the sphere at infin- 
ity. Note the apperance of (as opposed to G N ) in (3.1) - this ensures the correct weak 
field/slow- motion limit of the Einstein-asther theory, as we will see below. 

We can further express the right hand side of (3.1) in terms of the surface gravity on 
the sphere at infinity following (2.19). Using the asymptotic expressions of (4.6) - in which 
a particular choice of coordinates have been made - in (2.20), we find 



k/ dA( "- S,( *-* ) = 4^/ d ^- s) =2^- (32) 



M AD „ |r( . 

23oc 23 

Our starting point of the derivation of the Smarr relation is equation (2.30). As already 
noted, structurally (2.30) resembles the source free Maxwell's equations with T a b akin to a 
purely electrostatic field. In particular, by Gauss' law the flux of J- ab through the sphere at 
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asymptotic infinity equals the flux through sphere at the universal horizon 14 . Performing the 
flux integrals and using the expression for the ADM mass (3.2), we arrive at the promised 
Smarr relation 

where A im is the area of the universal horizon, q VH is the value of q (2.31) on the universal 
horizon, 

q VH = (1 - ci 3 )k uh + ^-K vn \x\im , (3.4) 

and we have used the expression (2.21) for the surface gravity at the universal horizon, k uh , 
to write q VH as above. Following [31, 32] and [18] we can furthermore introduce given by 

M 8 =(l- °-f) M ADM , (3.5) 

which is the total mass of an asymptotically flat solution defined in the asymptotic aether rest 
frame. Using (3.2) and the relation (2.6) between G N and Gas, we then have [18, 31, 32] 

Mae = Mc^n = M AOM G^ . (3.6) 

ZLx N 

The above relation between and ro ensures that one gets the correct Newtonian limit of 
the Einstein-aether theory far away from the sources [26, 31, 32]. In terms of the total mass, 
one also obtains a more natural presentation of the Smarr formula (3.3), namely 

iTrG^ 

To obtain the first law for the Einstein- aether black holes, we need to consider a variation 
which takes us from a given regular solution to a distinct nearby regular solution. The 
key observation leading to the first law is that the regular Einstein-asther black hole solutions 
depend on the single dimensionful parameter ro introduced in (2.32). As a result, the location 
of the universal horizon, r UH , is related to ro through a relation of the form r UH = ^ro, where 
/x is a dimensionless quantity which can depend only on the coefficients C2, C13 and C14. 
From the proportionality between r UH and ro, we now have q VH ~ 1 from (2.32), while 
A VH = 4-7rru H ~ and hence 5q im A VH = — ^uh^Ajh- Considering then a variation of the 
Smarr relation (3.7), the first law for Einstein-aether black holes follows in a straightforward 
manner [33] 

= . (3.8) 

Note that our derivation of the Smarr formula and the first law makes it manifest that 
at least when spherical symmetry is present, we can always have a "first law" applied to any 



14 According to (A. 10), the flux of Tab through any two-sphere, 23 r , at radius r, is the surface integral of q 
over 23 r . By spherical symmetry, this flux is the value of q at r, i.e., q(r), times the area of the sphere S,. 
itself. 
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sphere , where a variation of the total mass is proportional to q evaluated on the sphere, times 
the variation of the area of the sphere. Furthermore, since on dimensional grounds q ~ r^ 1 
as well asK~ r^ 1 , where q and k are the values of the respective quantities evaluated on the 
sphere in question, we can always write such a first law in terms of the surface gravity on the 
sphere. However, the importance of the universal horizon rests on the fact that it is the causal 
boundary in the spacetime and therefore, only (3.8) should possibly have a thermodynamic 
interpretation. 

4 The solutions 

In this section, our main goal is to present two exact, asymptotically flat, static, spherically 
symmetric, single parameter families of aether black hole solutions. These solutions provide 
additional evidence for the general result [21, 22] that all asymptotically flat, static and spher- 
ically symmetric aether black holes depend on a single parameter after imposing a regularity 
condition (to be discussed below). This single-parameter dependence, as already emphasized 
earlier, is crucial in our derivation of the first law. With our exact solutions, we furthermore 
verify the Smarr formula (3.7) and the first law (3.8). As we will also see, the interesting piece 
in q UH (3.4), that depends on the trace of the extrinsic curvature at the universal horizon, is 
absent (for separate reasons) for both these special solutions. 

In the following, we first adopt a convenient coordinate system to express the equations 
of motion (see the paragraph following (2.29)). Next, we present an asymptotic series solution 
of these equations, valid for large r and for arbitrary nonzero values of the couplings eg, C13, 
C123 and C14. The asymptotic solution has already been obtained in previus work [21, 22], and 
our purpose of presenting it here is three-fold: First of all, the asymptotic analysis determines 
the asymptotic (and sometimes the exact) nature of various relevant functions in the problem 
(e.g. the functional form q in (2.32)) and allows us to obtain the ADM mass (3.2). Secondly, 
the asymptotic analysis reveals that the general aether black hole solution can depend on 
at most two parameters, thereby providing a natural route to the topic of regularity of the 
solutions. Finally, the asymptotic analysis for two special choices of coupling constants lead 
to the exact solutions mentioned above. We discuss these special solutions in subsections 4.1 
and 4.2, respectively. 

To perform an asymptotic analysis of the equations, we set up an Eddington-Finklestein- 
like coordinate system which naturally respects the symmetries of the problem. With this 
choice of coordinates, the metric takes the form 

ds 2 = -e(r)dv 2 + 2f(r)dvdr + r 2 d^ , (4.1) 

and the timelike Killing vector is 

X a = {1,0,0,0} . (4.2) 

15 This point has also been stressed in [21], where a first law for the aether black holes applied to the spm-0 
horizon was obtained. For earlier work on the first law for an aether black hole at the Killing horizon using 
the Noether approach [34], see [32]. 
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The asther field can be parametrized as 

u* = {a(r),/?(r),0,0}, (3(r) = e( ffi^ \~ 1 , (4.3) 

2f(r)a(r) 

where the relation between a(r) and /3(r) takes care of the unit norm constraint (2.4). There- 
fore, to perform the asymptotic analysis, we only need the asymptotic behaviour of the three 
functions e(r), /(r) and a(r), which are given as follows 

e(r) = l + O^ 1 ), f(r) = l + 0{r~ 1 ), a(r) = 1 + 0(r _1 ) ^ (3{r) = 0(r -1 ) . (4.4) 

The boundary conditions on the metric coefficients are such that the solution is asymptotically 
flat, while those for the aether components are such that 

lim u a = {1,0,0,0} . (4.5) 

r— voo 

It is worthwhile to make the following comments at this stage: we set up the asymptotic 
analysis for arbitrary non-zero values of the coefficients C2, C13 and C14 such that C123 ^ 0. 
Although the results we present do not show this explicitly, if one takes the limit C123 — > 
of this solution, most of the coefficients in the 1/r expansions of e(r), /(r) and a(r) diverge. 
Therefore, we need to set up a separate asymptotic analysis when C123 = (but C14 / 0), and 
we then discover the exact solution presented in section 4.2. On the other hand, if one takes 
the C14 —7- of the general asymptotic solution by keeping C123 non-zero, the exact solution 
presented in section 4.1 is obtained 16 . 

We can now solve the Einstein's equations and the aether's equation of motion order 
by order in 1/r. We have solved the relevant equations to 0(r~ 10 ) which determines the 
functions e(r), f(r) and a(r) to 0(r~ 8 ) 17 . Not surprisingly, the asymptotic forms of these 
functions, as well as those which depend on them are quite cumbersome and do not convey 
too much information beyond that they can be found. We thus quote the relevant results 
only up to 0(r~ 2 ). To begin with, the metric components are 

e(r) = 1 " 7 + 0(r~ 3 ), f(r) = 1 + ^0 + 0(r~ 3 ) , (4.6a) 
while the components of the aether are 

a(r) = 1 + g + + 0(r-% /3(r) = -0 + 0{r^) . (4.6b) 

Our results agree perfectly with those in [22] (see equations 24—26), under the identification 
^1 «-> -r and A 2 «->■ (frg - r|). 

16 The scalar equation of motion (2.29) becomes trivially satisfied if both C123 and c\a vanish simultaneously, 
and consequently the structure and solution space of the field equations changes significantly. 

17 A11 the equations are second order ODEs in r and hence at 0(r~' n+2 - ) ) the functions only up to 0(r~ n ) 
contribute. 
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Given the results in (4.6a) and (4.6b), we can compute the series expansion for everything 
else; for example the asymptotic forms of (a ■ s), (u ■ x) an d (s • x) are 



(as) = ^+0(r- 3 ), (u-x) = -l+^ + ^+0( r - 3 ), (s- X ) = ^+0^) , (4.6c) 
Zr 2r or r 

respectively. The various components of the extrinsic curvature K a ^, as well as its trace, are 
likewise 

K Q = 2 ^ + 0(r- 5 ), K = - 2 ^ + 0(r~ 5 ), K = 0(r- 5 ). (4.6d) 

These results are useful to compute the ADM mass (3.2). 

The solution, at this stage, depends on two parameters namely ro and r^. Among these, 
the length scale ro is akin to the "Schwarzschild radius" and is related to the total mass of 
the black hole according to (3.6). The parameter r m is essentially the 0(r~ 2 ) coefficient of 
a(r), and is defined in this way for convenience. From the asymptotic analysis, it appears 
that is a second free parameter on which the solutions depend. This, as we now explain 
following [21, 22], is not the case after all; rather r^ is related to ro upon requiring that the 
solutions are regular everywhere outside the universal horizon. 

We begin our discussion with the observation that the Einstein-aether theory admits a 
"ground state" solution where the spacetime is four dimensional Minkowski and the asther is 
{1,0,0,0} with respect to an observer in the preferred frame (called the asther rest frame). 
In [15], the authors consider perturbations around this background and show, in particular, 
that there is a spin-0 18 mode which propogates with a speed sq given by 

4 = n Cl23 t~ CU) , o s ■ ( 4 - 7 ) 



cu(l - ci 3 )(2 + ci 3 + 3c 2 ) 



with respect to the asther rest frame. Because of general covariance of the Einstein- asther 
theory, perturbations around an aether black hole background will also give rise to a spin-0 
mode with a local speed given by (4.7). The spin-0 horizon is a hypersurface beyond which 
any outward moving excitation travelling with sq (or less) gets trapped. More precisely, the 
spin-0 horizon is hypersurface where the timelike Killing vector becomes null with respect 
to the "effective spin-0 metric" gj^ = g a b — (sq — l)u a Ub [21, 22]. Equivalently, we can also 
define the spin-0 horizon as the hypersurface where (s • x) 2 = Sq(u ■ x) 2 - 

For generic values of the couplings C2, C13 and C14, which respect (2.5), the spin-0 speed 
so is a non-zero finite quantity. Consequently, the spin-0 horizon can be located anywhere 
outside the universal horizon. However, for finite non-zero sq one can always apply the field 
redefinitions introduced in [35] to set so = 1, thereby making the spin-0 horizon coincide with 
the Killing horizon [21, 22]. This extra condition therefore reduces the number of independent 
couplings from three to two. However, we should emphasize that this does not mean we are 

18 The results of [15] show that when the aether is not hypersurface orthogonal and if no symmetry is assumed, 
there are additional spin-1 and spin-2 modes. In our case, the condition of hypersurface orthogonality of the 
aether will prevent any spin-1 mode from propagating in the black hole backgrounds we consider. Likewise, 
any spin-2 mode will be excluded because of spherical symmetry. 
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exploring a restricted coupling space. Rather, we are using an extra freedom in the theory 
(the field redefinitions) to conveniently choose (by imposing sq = 1) typical sets of couplings 
{c2,ci3,ci4} which label larger equivalent classes. In [21, 22], the authors use the above logic 
to effectively scan a smaller coupling space in their numerical constructions of asymptotically 
flat, static and spherically symmetric aether black holes. Those studies clearly prove that 
a solution for generic values of tq and r m is singular, precisely at the location of the spin-0 
horizon. However, once the solution is required to be regular everywhere outside the universal 
horizon, the extra constraint automatically makes r m dependent on Tq, i.e., the former cannot 
be an extra parameter. Since the general asymptotically flat solution can at most depend 
on two parameters, the regularity condition reduce the parameter space so that we have a 
one-parameter family of solutions. In this manner, [21, 22] obtain a unique asymptotically flat 
asther black hole solution for a given value of the parameter rn (and a given set of couplings) 
by making the solution regular at the corresponding spin-0 horizon. 

The field redefinitions of [35] (see also [21]) mentioned above also scale C123 and (1 — C13) 
in the same way while keep C14 invariant. It is then clear that such a transformation does not 
exist when either of C123 or C14 vanish or when C13 = 1 (we can rule out this last possibility 
owing to (2.5)). At the same time, according to (4.7), the spin-0 speed diverges as C14 — > 0, 
while it vanishes as C123 — > 19 . In the context of a black hole solution, when C14 = 0, the spin- 
horizon coincides with the universal horizon, since as noted in the introduction, the latter 
is the causal boundary for arbitrarily fast excitations. On the other hand, when C123 = 0, the 
spin-0 horizon in a black hole solution is pushed all the way to spatial infinity and so overlaps 
with the asymptotic boundary. Owing to the absence of the field redefinitions for these cases 
however, the spin-0 horizons cannot be mapped on to the metric horizon. Remarkably, there 
exists exact solutions for these special cases, where the spin-0 regularity condition can be 
illustrated in an explicit manner. We present these solutions in the following two subsections. 

4.1 Exact solution for C14 = 

When the coupling C14 is set to zero, the system admits an exact solution given by 

e(r) = 1 " 7 " f(r) = 1 , (4.8a) 

and 

«W = ^y + * M = -%- Mb) 

From the explicit solution, we further work out 

(s-x) = -P(r) = r fi, (" • X) = ~Ve(r) + W? = -\A ~ 7 + ^ ~ ■ ( 48c ) 

19 We note that there are other limits of (4.7) when so can vanish or diverge: C14 — ¥ 2 (so vanishes), C13 — ¥ 1 
(so diverges) and (2 + C13 + 3C2) — ¥ (so diverges). However, they all violate the constraints (2.5), and therefore 
are excluded on physical grounds. 
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As mentioned in the beginning (2.5), we always assume 5= C13 < 1. 

We now investigate the locations of the Killing and universal horizons in this solution. By 
definition, the Killing horizon is located where (x'X) = 0> or equivalently, at the largest root 
of e(r) = 0, and the universal horizon is located at the largest root of (u-x) = 0. From (4.8a) 
and (4.8c) this amounts to solving two quartic equations. Rather than doing this directly, we 
extract most of the important properties of these roots from simple arguments. 

We begin by noting that e'(r) > everywhere owing to C13 > 0. Therefore, e(r) is a 
monotonically increasing function and we have a single real root 20 at r = r KH , which is the 
location of the Killing horizon. Secondly, from (4.8c) (u • x) 2 = e ( r ) + ( s ■ x) 2 ~ therefore 
from the monotonicity of e(r), even the largest root of (u-x) 2 is necessarily located at some 
T = t*uh < r i<H- We also conclude that (u-x) 2 = has at most two real roots by noting that 
the function has a single minimum at r = ^/4(1 — c\^)r^/rQ. Furthermore, the two roots 
are distinct when < r* , they coincide when = r^, and there are no real roots when 
fst>f*, where 



4 



27 



1/4 



.1 ~ c 13. 

The situation should be contrasted with the existence of an event horizon for the usual 
charged (Reissner-Nordstrom) black hole. However, there is an important difference between 
the present solution and the charged black hole solution in regards with the regularity of the 
solutions everywhere. In case of a charged black hole, the solutions are regular everywhere 
except at r = and are physically allowed as long as the extremality condition is met. To 
examine the regularity of the present solution, we can see from the expressions of the curvature 
scalars for the present solution 

that the ambient spacetime is free of any curvature singularities except at r = 0. But, 
the aether field being a physical component of the theory, we also need to make sure that 
the solution for the aether is regular everywhere as well. A coordinate independent quantity 
associated with the aether, which can signal the existence of pathologies in the present solution 
is (u - x) 21 ■ Indeed, when r < r UH , (u ■ x) 2 is negative based on our discussion above, and 
hence (u ■ x) is purely imaginary between r = r UH (the location of the universal horizon) and 
the smaller root of (u ■ x) 2 = 0. In other words, for generic values of r^, the aether solution is 
irregular at the universal horizon. Naturally, this irregularity is prevented if (u-x) 2 is never 
allowed to be negative. This regularity condition, along with the demand for the existence of 
at least one root 22 of (u ■ x) 2 = 0, uniquely implies that the regular physical solution exists, 

20 Of course, e(r) = must have at least two real roots, but the second real root must be negative by the 
above argument, and hence unphysical. 

21 (s • x) f° r this solution is another coordinate independent quantity. But it is nicely behaved everwhere 
except at r — 0, and is therefore incapable of signaling irregularities. 

22 As required by cosmic censorship - since we have superluminal propagation, the Killing horizon cannot 
save cosmic censorship. 
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Figure 1. r KH in units of r as a function of C13. 
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nl/4 

(4.9) 

.1 - C13 

Thus, Tgg is not an independent parameter. We have already argued that the spin-0 horizon 
for this solution overlaps with the universal horizon. Therefore, the regularity condition is 
indeed a regularity condition at the spin-0 horizon. 

Manifestly, the regular solution depends on a single parameter ro. Here onwards, when 
we talk about this exact solution as well as about any quantity pertaining to it, the condition 
(4.9) will always be implied. The location of the universal horizon for this physical solution 
is very easy to find - it is a root of both (u • x) 2 = and d/dr(u • x) 2 = 0, and is given by 

r UH = ^ • (4.10) 

Quite interestingly, the result does not depend on the value of C13. The location of the 
Killing horizon r KH however does depend on C13, but we did not attempt to find an analytical 
expression for it; instead, we solved for r KH numerically and the result is presented in figure 1. 

Prom its definition (2.24) T^. vanishes when C14 = and the aether stress tensor becomes 
diagonal, with the nontrivial components given by 

3c 13 rj 3c 13 rj > 12c 13 rj 

From the equation of motion (2.24) we then have V S K = 0, i.e., K is constant on a given 
hypersurface £[/. But K vanishes asymptotically for asymptotically flat spacetimes - this can 
be most readily seen from the fact that u a ~ \ a asymptotically, so that K ~ V • x = 0. Thus 
K = on every hypersurface "Ejj, an d therefore everywhere in spacetime. A related point is 
that, manifestly, the solution does not depend on the coupling c% in any way. According to 
[28], C2 is the coupling of the K 2 term of the aether Lagrangian. Therefore we see that every 
reference to C2 has been removed as K vanishes. 

We are now in a position to derive a Smarr relation. Given the solution (4.8) we can 
compute the surface gravity on the universal horizon following (2.21) and this turns out to 
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be 

^uh = ^ 7T 8 r • (4.12) 

9r (l - ci3) 

Using A VH = 47rru H = 97rrQ/4 as the area of the universal horizon, we can therefore derive a 
Smarr formula for the present solution 23 



(1 - ci 3 )k uh ^ 



(4.13) 



Varying the parameter ro we then also obtain a first law of black mechanics for the present 
solution explicitly 

(1 - ci 3 )k uh SA vh 



8vrG a 



(4.14) 



Comparing (4.13) and (4.14) with the general Smarr relation (3.7) and the first law (3.8) re- 
spectively, we find perfect agreement. Naturally, the interesting piece proportional to ci23-Kuh 
is absent as K = everywhere in this solution. 

4.2 Exact solution for C123 = 

When the coupling C123 is set to zero, we also have an exact solution given by 



e(r) = l- — ' /( r ) = 1 > «( r ) 

where r u is a positive constant, given by 



1 + 



2r 



cu 



2(1 -ci 3 ) 



2 



(4.15) 



(4.16) 



The requirement of positivity of r u follows from demanding the function a(r) be regular 
everywhere. Consequently, this imposes the following bound on the couplings C13 and C14 



C14 S 2ci 3 • 



(4.17) 



This is consistent with the constraints ^ C13 < 1 and ^ C14 < 2 (2.5), as we assume in 
general. The positivity of r u criterion also rules out another possible solution for r u 



2(1 - ci 3< 



+ 1 



ro 
2 ' 



which solves all the equations of motion, but is manifestly negative for all values of C13 and 
C14. The curvature invariants for this solution are 



m = 0, 



ab 



4r 2 (r + r u ) 



Note: M m = M ADM for this solution since C14 = 0. 
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demonstrating that the geometry is free from any curvature singularity everywhere except at 
r = 0. Thus the solution is regular everywhere. 

The reader might have spotted that the parameter r m does not appear in this solution. 
Instead we have the parameter r u , which is however not a free parameter. Rather, there 
are two possible choices for r u as functions ro and this ambiguity is resolved (in the form of 
choosing r u positive) by demanding that the solution be regular everwhere except at r = 0. 
One way to appreciate the difference between the present solution and all the solutions with 
C123 7^ 0, is to note that we need separate asymptotic analysis for the cases C123 7^ and 
C123 = 24 . As it turns out, when C123 7^ 0, the equations of motion force the C(r _1 ) 
coefficient of a(r) to be ro/2 (4.6b). When C123 = 0, this requirement no longer holds and we 
are left with a free parameter r u . Furthermore, when C123 7^ 0, the parameter r m appears as 
the 0(r~ 2 ) coefficient of a(r) and stays as the second free parameter in the general asymptotic 
analysis. For the C123 = case however, the 0(r~ 2 ) coefficient is a function of r u , and as 
the subsequent analysis reveals, the equations of motion eventually restrict r u to take one of 
the two choices mentioned above. The requirement for regularity then removes one of the 
choices. We have already argued earlier that the spin-0 horizon for this solution "coincides" 
with the asymptotic boundary at infinity. We thus see that even with the correct boundary 
conditions imposed, there can be two different solutions corresponding to two different values 
of the paremeter r u 25 . The actual regularity condition here comes in the form of choosing 
the correct value of r u as given in (4.16). Manifestly, the regular solution depends on a single 
parameter ro. 

We next adress the issue of the locations of the Killing and universal horizons. From the 
explicit solution, we further work out 

{s ■ x) = -)8(r) = („. X ) = -! + £. (4.18) 

The roots of e(r) = can be explicitly found in this case and they are r = (ro + r u ) and 
r = — r u respectively. The second root is negative, i.e. unphysical, and the unique Killing 
horizon is at 

r K H = r + r u . (4.19) 

On the other hand, (u ■ x) has a unique root at r = ro/2 which is therefore the location of 
the universal horizon in this solution 

r UH = V j . (4.20) 

As in the case of the exact solution with C14 = 0, the location of the universal horizon does 
not depend on the coupling constants for the present case either, while the location of the 
Killing horizon does. 



24 There is no such need when the coupling in question is C14 as the general asymptotic solution (4.6) admits 
a smooth C14 — > limit. 

25 The choice of the correct asymptotic boundary conditions with appropriate fall off conditions for the 
various fields should be treated as partly making the spin-0 horizon regular in this case. In particular, this 
seems to be the reason for no additional continuous parameter in the solution. 
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With vanishing according to (2.29) the stress tensor is also diagonal in this solution, 
with the non-trivial components given by 

rr 33 = 2c ^(ro + 2r u ) 2 -c M r 2 ~ 

To derive a Smarr relation, we proceed as before and compute the surface gravity on the 
universal horizon from (2.21) using the present solution (4.15) 

« UH = n 2 ~ C1 ' • (4.22) 
(1 - ci 3 )r 

This time, using A VH = 47rry H = nr^ as the area of the universal horizon, we arrive at the 
Smarr formula for the present solution 

Mae = (1 " C ; 3 ^ UH ^ UH . (4.23) 

Therefore, the first law of black mechanics for the present solution is 

SM m = 

We find perfect agreement, once again, upon comparing (4.23) and (4.24) with the general 
Smarr formula (3.7) and the first law (3.8) respectively. This time, the interesting piece 
proportional to c1.23.KuH is absent for the obvious reason. 



5 Summary and Conclusions 

In this paper we have studied static, spherically symmetric, asymptotically flat black hole 
solutions of the Einstein-aether theory, a generally covariant modification of general relativity 
where a vector field, the aether, is forced to satisfy a unit normalization constraint. Since 
there is a preferred frame of reference defined by the aether, the solutions of the theory 
violate local Lorentz invariance and so matter fields do not necessarily have a finite local 
limiting speed. Even though at first sight the notion of a black hole seems impossible in 
such a situation, earlier work [21, 22] provided concrete evidence in support of the existence 
of a single parameter family of static, spherically symmetric, asymptotically flat black hole 
solutions. The causal boundary in question, called the universal horizon, is a hypersurface 
which traps excitations travelling even at infinite local speed. 

In this work, we extend these earlier studies by demonstrating that a Smarr relation and 
a first law of black hole mechanics associated with the universal horizon can be found. We 
also provide analytical evidence for the existence of a universal horizon by constructing two 
exact solutions for special choices of the couplings of the Einstein-aether theory, and verify 
the Smarr formula and the first law with these special cases. Critical to our proof of the 
first law is the fact that the solutions depend on a single parameter, namely the mass of the 
solution [21, 22]. 
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The Smarr formula (3.7) and the first law (3.8) suggest that in the present context, the 
quantity q VH (3.4) plays the same role as surface gravity at the Killing horizon of a black 
hole in general relativity. Indeed, from (3.4) q VH does include a contribution from the surface 
gravity at the universal horizon, but there is also an additional piece that depends on the 
extrinsic curvature of the preferred foliation as it asymptotes to the universal horizon. Based 
on the causal nature of the universal horizon, a thermodynamic interpretation of the first law 
(3.8) seems necessary, along the lines of [2]. However, a complication arises as in Einstein- 
aether theory Cerenkov radiation generically occurs since aether-metric modes and matter 
fields all have different speeds. In particular, a matter field mode propagating outwards from 
the universal horizon would emit spin-0 Cerenkov radiation, thereby modifying any thermal 
spectrum. Since no detailed examination of the radiation spectrum from a universal horizon 
has been made yet, it is unclear whether or not there is a Cerenkov-type component in addition 
to the usual thermal flux. Note, however, that when cu — > the speed of the spin-0 mode goes 
to infinity, and so there is no possibility of emission of spin-0 Cerenkov radiation. On the other 
hand, when C123 —> 0, the speed of the spin-0 mode goes to zero, and so any spin-0 Cerenkov 
radiation would carry no energy away from the universal horizon. Thus, in both cases the 
reduction to a situation where the presumed temperature is proportional to the surface gravity 
is consistent with a limit where any Cerenkov radiation would naturally become unimportant. 
For now, we leave the detailed examination of whether or not a thermodynamic interpretation 
of q ull exists to future work. 

Acknowledgments 

We are grateful to Ted Jacobson for critically reviewing an earlier draft of this paper and for 
providing numerous invaluable suggestions. P.B. acknowledges the hospitality of the KITPC, 
Beijing, the Berkeley Center for Theoretical Physics and the theory group at CERN. The 
work of P.B. and J.B. is supported by the NSF CAREER grant PHY-0645686. D.M. thanks 
the University of New Hampshire for research support. 

A Some algebraic details 

In this appendix, we provide some of the technical algebraic details behind our analysis. Of 
central importance to our analysis is the hypersurface orthogonality relation (2.11), as well 
as the assumptions of spherical symmetry and staticity. We already noted in section 2 (see 
the footnote around (2.17)) that as a consequence of the hypersurface orthogonality of u a , 
spherical symmetry and staticity, the vector s a is hypersurface orthogonal with respect to 
the foliations S s of the spacetime, and satisfy s A ds = 0. Contracting these hypersurface 
orthogonality conditions with u a and s a respectively, we obtain 

V[ a n b] = -(a • s)u [a s b] , V [a s b] = -K u [a s b] (A.l) 

Further contractions of these relations with any Killing vector rf and with u a and/or s a 
lead to extremely useful scalar identities which play a major role in our analysis, and in 
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particular, in the derivation of (2.30). One of our goals in this appendix is to briefly outline 
how these identities are obtained in a generalized and coherent fashion. To that end, we recall 
(as already noted in the paragragh following (2.15)) that by the spherical symmetry of the 
problem, any rank-two tensor can have components only along the bi- vectors u a u b , ur a Sj,\, 
u [a s b]i s aSb and g a b. The most basic among these basis-expansions are those for V a Ub 

k 

V a ii b = -(a- s)u a s b + K ab , K ab = K s a s b + —g ab (A. 2) 

and V a Sb 

V a s b = K s a u b + K$, K$ = -(a- s)u a u b + ^g ab (A. 3) 

(s) 

where K\' is the extrinsic curvature of the hypersurfaces S s due to their embedding in the 

spacetime, and k and K are the traces of the extrinsic curvatures of the two-spheres 23 due 
to their embeddings in T,jj and T, s respectively 

k = (l/2)g ab £ sgab , K = (l/2)g ab £ ugafc (A.4) 

Our results in this paper make heavy uses of (A. 2) and (A. 3). 
Now, consider an arbitrary vector of the following form 

A a = -fu a + hs a (A.5) 

where / and h are arbitrary functions respecting the symmetries of the problem. A natural 
construct is the two-form 

Fab = V [a A b] = Q u [a s b ] (A. 6) 

where the second equality follows by the spherical symmetry of the problem, and the scalar 
Q is given by 

Q = f(a ■ s) + V,/ - hK - V u h (A.7) 

Using the torsion-free condition of the covariant derivative, it can be shown that Q, like F ab , 
is invariant under the "gauge transformations" A a ^ A' a = A a + V a A. 

Comparing (A.l) with (A. 6), we see that the former relations are special cases of the 
latter. Contracting (A. 6) with a Killing vector rf yields 

V a (A-r ) ) = Q(s-r})u a -Q(u-r t )s a (A.8) 

Contracting (A.8) further with u a and s a gives 

V u (A- V ) = -Q(s-rj), V s (A-rj) = -Q(u- V ) (A.9) 

We make ample use of the special cases of (A.8) and (A.9) for A a = u a and A a = s a 
throughout the analysis. 
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The flux of F a b over any two-sphere 23 is 



/ 

s 



dZ ab F ab 



f dAQ 

B 



(A.10) 



Since F ab is antisymmetric, we have a kinematically conserved current J a through 



V 6 F 



«6 



(A.ll) 



By spherical symmetry, the current can only have components along u a and s a . These pro- 
jections are given by 



u ■ J = -(kQ + V S Q) = -V • [Qs\, s-J = -(KQ + V U Q) = -V • [Qu] 



(A.12) 



where v 0, = F a b v b and v a = F a b v b are projections of any four vector v a on the hypersurfaces 
S{7 and S s respectively, P a fe = u a Ub + 5 a b and P a 6 = — s a Sb + 5 a b are the corresponding 
projection tensors, and V a = IP^Vb and V a = !P a b Vb are the corresponding projections of 
the covariant derivatives. 

The most important example of a kinematically conserved current, which is central in 
the derivation of (2.30), is that due to the timelike Killing vector x a - Choosing the vector 
A a = (l/2)xa in (A. 5), so that from (A. 6) F a b = V a Xb, we obtain Q = —k. From the identity 
VfcV a x b = 9L a bX b i which forms the basis of the Smarr relation, we find that the kinematically 
conserved current is 



J^ = 3l a b X , u- J X = V -{ks\, s-J x = V-[ku] 
We conclude this appendix by noting the following two identities 



(u ■ rf) (V • v) = V ■ I — (u ■ v)(s ■ 77) + (s • v)(u ■ rj) |i 



and 



(s ■ rj)(V • v) 



(u ■ v)(s ■ rf) + (s ■ v)(u ■ r?) j 



u 



(A.13) 



(A.14) 



(A.15) 



valid for any four-vector v a . The identities (A.14) and (A.15) are indispensible for dealing 
with the contributions from the asther in the derivation of (2.30). 



B More on spin-0 horizon (ir) regularity 

In this appendix, we present some further studies of the equations of motion, where the 
irregularity of a general asther black hole solution near its spin-0 horizon for generic values of 
ro and [21, 22] manifests itself in interesting ways. 

We have performed a near-universal-horizon analysis of the equations of motion to lowest 
order in the near-universal-horizon radial coordinate. The most striking feature of the analysis 

26 We call J a a kinematically conserved current if it satisfies the continuity equation V • J = irrespective 
of any equations of motion. 
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is that two separate treatments are required for the cases C14 7^ and C14 = 0, since, the series 
coefficients are divergent in C14 in the C14 7^ case 27 . Indeed, in this case the near-univeral- 
horizon region is also the near-spin-O-horizon region when C14 is sufficiently small. Since we 
did not impose any a priori regularity condition on the solution in our analysis, this divergent 
behaviour is clearly due to the irregularity of the solution at the spin-0 horizon. 

The near-universal-horizon analysis for C14 = 0, on the other hand, successfully repro- 
duces 28 the near- universal-horizon behaviour of the exact solution for C14 = 0. Likewise, upon 
taking the limit C123 — > of the near-universal-horizon solution for C14 7^ 0, we successfully 
reproduce the near-universal-horizon behaviour of the C123 = exact solution. 

We have also studied perturbative corrections to the exact solutions for C14 = (4.8) and 
C123 = (4.15) - for small C14 in the former case, and for small C123 in the latter. The former 
analysis (i.e around the C14 = exact solution) did not yield a complete solution because 
of the very complicated nature of the relevant equations - we could only find the £)(c\±) 
correction to the function f(r). In the context of the C123 = exact solution however, with an 
additional mild assumption, namely cu = 2ci3 — e where C123 = e < 1, we found the 0(0123) 
correction to all the three basic functions e(r), /(r) and a(r) 29 . The most important property 
of the perturbative corrections, relevant for the present discussion, is that the corrections have 
logarithmic singularities at the location of the unperturbed spin-0 horizon 30 . More specifically, 
the 0(0123) correction functions to the C123 = exact solution depend on log(r/ro) 31 , while 

27 The situation should be compared with the asymptotic analysis when C123 7^ and the corresponding 
divergence of the asymptotic solution in C123 - see the comments following equation (4.5). As discussed earlier, 
for small C123 the spin-0 horizon is pushed towards asymptotic infinity. Therefore the asymptotic solution is a 
good approximation to the near-spin-O-horizon region of the solution for sufficiently small values of C123. The 
divergence here is therefore a manifestation of the irregularity of the general asymptotic solution (where no 
regularity is imposed) at its spin-0 horizon. 

28 Of course, the exact values of the functions at the universal horizon, as well as the radial location of the 
universal horizon in the Eddington-Finklestein radial coordinate (4.1), depend on the normalization of the 
asymptotic data and the actual definition of the radial coordinate itself, and cannot be "derived" in this kind 
of an analysis. 

29 In our analysis, the spin-0 speed so can be kept arbitrarily small as long as C123 <C C13 < 1; therefore, 
the spin-0 horizon can be "arbitrarily far away", i.e., "arbitrarily close" to the asymptotic boundary. Looking 
at the expression for r u (4.16) we further notice that the assumption on C14 implies that the perturbative 
correction we are studying can also be thought of as a perturbation around a Schwarzschild background. We 
also refer the reader to a useful study of non-spherical perturbations about black hole solutions in extended 
Hofava gravity and Einstein-asther theory [36]. 

30 Even though our analysis of the small C14 perturbation around the C14 = exact solution is incomplete, 
the statement is true for the &(ci4) correction for the function f(r) - see below. The relevant equations in 
this context also indicate that such logarithmic singularities will be present for e(r) and a(r). 

31 One might be worried that the 0(0123) correction functions to the C123 = solution destroy the asymptotic 
flatness condition, due to the presence of the log(r/ro) in their expressions. We however do not claim that the 
C(ci23) correction approximates the actual solution everywhere in spacetime. Indeed, the 0(0123) correction 
can be a good approximation near the universal and the Killing horizons, but we need to consider higher order 
corrections in C123 to approximate the actual solution as we move out in r, and presumably we need to go to 
all orders in C123 in order to "approximate" the actual solution at asymptotic infinity. On the other hand, by 
the results of [21, 22] we are guaranteed to reach an asymptotically flat regular solution (provided we choose 
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for the C14 = exact solution, the 0{c\/^) correction to f(r) depends on log(l — 3ro/4r). Now, 
although the exact solutions are regular everywhere, their first order corrected counterparts 
are a priori not so. The logarithmic divergence at the unperturbed spin-0 horizon therefore 
clearly signals an irregular behaviour 32 . 
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